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We discuss a general form of the exchange energy for a homogeneous system of interacting elec- 
trons in two spatial dimensions which is particularly suited in the presence of a generic spin-orbit 
interaction. The theory is best formulated in terms of a generalized fractional electronic polariza- 
tion. Remarkably we find that a net generalized polarization does not necessarily translate into an 
increase in the magnitude of the exchange energy, a fact that in turn favors unpolarized states. Our 
results account qualitatively for the findings of recent experimental investigations. 

PACS numbers: 71.10.Ca, 71.10.-w, 71.70.Ej, 73.61. Ey 
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In most cases, as a first step in assessing the role of the 
electronic interactions, an analysis of the role played by 
the exchange energy is illuminating and most often very 
useful. As we show here, the exchange energy acquires 
particularly interesting and unusual properties in two di- 
mensions in the presence of spin-orbit of the Rashba type. 
In this case the chiral nature of the single particle states 
involved competes with standard requirements that are 
commonly identified with the exchange interaction. 

The interplay of spin-orbit coupling and Coulomb in- 
teraction is unexplored, except for the initial and limited 
work of Refs. [HQ where some quasi-particle properties 
were investigated. Most of the work is relevant in the 
high density limit in which only weak corrections to the 
standard behavior were identified. Systematic studies of 
the mean-field phase diagram of this interesting system 
are provided elsewhere 

Alongside the general theoretical relevance of such sys- 
tematic studies, a timely practical application of such 
an analysis can be found in the study of spin polariza- 
tion states in two dimensional hole systems in a GaAs 
heterostructure^ 

The expression for the exchange energy per particle in 
a homogeneous electron liquid is well knownji In partic- 
ular, for the two dimensional case, which is of relevance 
in the present discussion, the formula is given by: 



4\/2(l+p)l + 
37r r s 



(1) 



where we measure energy in Rydbergs and the signifi- 
cance of the sub index will be presently made clear. 
In Eq. (JTJ the parameter p is defined as follows: 



P : 



(2) 



where n is the particle number density and n± refers to 
the particle density corresponding to the two spin sub- 
bands. r s = \j y/ixa\n is the usual density parame- 
ter in two dimensions. It is important to notice that 
Eq. (JTJ is obtained under the assumption that the spins 
are quantized along a common arbitrary quantization 
axis z. Therefore, when, as in this particular case, the 



label ± refers to the orientation fl of the spin, p assumes 
the meaning of a fractional spin polarization. It is easy 
to see that the exchange energy (JTJ) monotonically at- 
tains its maximum magnitude for p = 1, and therefore 
favors the polarization of the system. In fact, within the 
Hartree-Fock approximation, a transition to a fully spin 
polarized state occurs for r s = . 37r — ~ 2.01.— 
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As we will presently show, this scenario changes in in- 
teresting ways in the presence of the spin-orbit interac- 
tion. In the quasi two dimensional electronic systems 
present in a high quality semiconducting heterostruc- 
ture, the latter can come in different guises ranging from 
Rashba and Dresselhaus couplings to terms induced by 
the application of an external magnetic fields The fol- 
lowing is a systematic study of the effects of spin-orbit 
interactions of various type on the exchange energy and 
the polarization properties of a clean two dimensional 
system. The theory can be elegantly formulated in a 
compact way by introducing a generic form of spin orbit 
coupling, something that allows us to draw quite general 
conclusions. 



The paper is organized as follows: Section [TJ introduces 
our model hamiltonian including a spin-orbit coupling 
of generic form; Section [TT| describes the solution of the 
corresponding non interacting problem in terms of the 
generalized polarization. In Section ITTT1 which represents 
the main part of the paper, the role of the exchange en- 
ergy is studied in detail. In particular, it is demonstrated 
there that the presence of spin-orbit coupling leads to a 
qualitatively different relation between this quantity and 
the generalized polarization. Section IIVI contains a dis- 
cussion of recent experiments on the polarization of holes 
in GaAs heterostructures in the light of these findings. 
Finally some general remarks and conclusions are pro- 
vided in Section fVl 
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I. MODEL HAMILTONIAN 



with energies^ 



We consider the following generic model two dimen- 
sional electronic system described by the hamiltonian 9 
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where the single-particle terms are of the form: 



(3) 



(4) 



the motion taking place in the x-y plane. In Eq. ((4]) we 
have defined p± = p x ± ip y and a± — a x ± ia y . Here n is 
an integer number assuming values from to 3. In view 
of its structure we will refer to the second term in ([4| as 
a generalized spin-orbit coupling. 

In the simplest case, n = corresponds to the familiar 
Zeeman coupling with 7 = ShsM. f or which the standard 
result (p} applies. More complicated n = expressions 
are contemplated in the Appendix. 

For n > 0, Eq. (H} describes different types of bona fide 
spin-orbit interactions. In particular, for n = 1 we obtain 
a form equivalent to the Rashb a 10 ! 11 or Dresselhausi^ 
spin-orbit hamiltonians, appropriate for two-dimensional 
conduction electrons in III-V semiconductors. Further- 
more, the cases n = 2 and n — 3 are appropriate for 
holes in III-V semiconductors as GaAs. Explicit expres- 
sions are provided in the Appendix. 

From a physical point of view the n — 2 term arises 
in the presence of an external in-plane magnetic field, 
the coefficient 7 being proportional to the value B of 
such field. For hole carriers in GaAs heterostructures, 
this term is the dominating low-field effect in the high 
symmetry growth directions [001] and [111], for which the 
Zeeman term can be safely assumed to be approximately 
vanishing 4 

Finally, the n = 3 term corresponds to a Rashba spin- 
orbit coupling for holes, the value of 7 being in this case 
approximately proportional to the value of the average 
electric field of the confining potential. 13 



£ k,± 



% z k 
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Notice that in view of the structure of the spinors the 
integer n can be seen as a spin direction winding number. 

It is not difficult to show 3 - that a class of determinantal 
many-body states with compact momentum space occu- 
pation which are homogeneous and isotropic in the plane 
of motion are uniquely determined by the areal density 
n and the fractional polarization p defined as in 
It must be immediately noted here that in the general 
case p should not be interpreted as a spin polarization 
for this is in general vanishing in these states, p does 
merely determines the two Fermi vectors k± of the so 
called spin-split subbands via the relations: 



k± = y / 2Trn(l±p) 



(8) 



which uniquely determine the occupation numbers rik± 
of the subbands. 

It is readily found that the non interacting energy 
per particle (in Rydbergs) of such Slater determinants 
is given by: 



°n V s,P) — o I .„ -, 1 n 



where we have defined the dimcnsionless coupling: 



m n-l e a(n-2) 

7= 7^-2 7 



(9) 



(10) 



The first term in (J9j) is the contribution of the kinetic 
energy while the second corresponds to the spin-orbit in- 
teraction. 

In this case, minimization of the total energy with re- 



spect to the value of p, 



at:'' 



0, leads to the (finite) 



ground state value of the equilibrium generalized polar- 
ization p^in ■ Interestinj 
dimensionless quantity 



ization p^ in ■ Interestingly the result only depends on the 



II. NON INTERACTING ELECTRONS 

The eigenstates of (|4|) are plane waves, a fact that 
allows one to write the spin-orbit term in the form 
7 (hk) n cr ■ Ski where is defined as follows: 



Sk = — sin(n0k)i + cos(n0k)y 



(5) 



Here </>k = arctan(fcj / /fc a; ) is the polar angle spanned by k. 
The unit vector Sk determines the direction of the quan- 
tization axis for the particular value of the wave vector 
k. The two possible spin orientations immediately give 
the eigenstates: 



</>k,±(r) 



±1 

p intf> k 



(6) 



g = 2 = 7^"" , 
so that it can be compactly expressed as follows: 
g for n = 



P 



(0) 



for n = 1 
for n = 2 

3 



for n 

For small coupling this quantity behaves as 



(0) 

Pmin 



n(n — 2) O r 
fj + ^ g V + 0(g 5 ) 



(11) 



(12) 



(13) 



One can then immediately notice that since for n = 1 
we have g = y2*fr s , in this case the high-density regime 



3 



is equivalent to a vanishing spin-orbit. The opposite ob- 
tains for n = 3, since g — 2v27/V s . For quadratic spin- 
orbit g — 27, and the fractional generalized polarization 
p is independent of the density. 

It is also useful to define here the depopulation cou- 
pling strength jjp as the particular value of 7 for which, 

at a given density, the upper band empties and p^in = 1 • 
This is readily obtained from (|12p to be given by: 



_(o) 
Id 
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for n = , 
for n = 1, 2, 3 



(14) 



occupation numbers nk± . The total exchange energy can 
be written in the following elegant general form: 3,4 



E r 
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El + n/j,' s k • Sk' 
Uk-k' ^ "k M n-kV , 

k,k / ;/x,/x / =± 

which represents a functional of rik± and Sk- Eq. (|15|) 
immediately reduces to the familiar textbook result when 
Sk = z. In this case only states with parallel spin con- 
tribute. 

Making use of Eq. in (fT5f leads to the following 
result for the exchange energy per particle (in Rydbergs): 



III. EFFECTS OF EXCHANGE 



We turn next to the effect of exchange. If translational 
invariance is not broken,— a single Slater determinant can 
be constructed with plane waves states characterized by 
generic orientation of the spin quantization axis Sk and 



£x,n(r a >p) = £xfl{r s ,p) + S£ x>n (r s ,p) , (16) 

where the correction to Eq. {1} can be obtained from the 
following quadrature: 



J 




x dx 
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V x2 + : 



2xy cos 9 



:A6 



(17) 



The resulting exchange energy is plotted in Figure [T] as 
a function of p. We notice that for p = 0, the exchange 
energy is independent of the spin quantization axis ori- 
entations Sk- This can be understood by realizing that 
the corresponding many-body state can be constructed 
by repeated application of 6 k+ 6 k _, an operator that cre- 
ates a spin singlet and is therefore independent of the 
spin quantization direction. As a consequence the s k de- 
pendence of all the physical quantities (e.g. the exchange 
and the spin orbit energies) stems only from the existence 
of regions of momentum space where nk+ 7^ ?^k- . It is 
also important to remark that only for n = is the mag- 
nitude of the exchange energy maximum for p = 1. For 
the other cases the minimum occurs at p — p\ ~ 0.915 
(although this is not obvious from Figure [T]) for n = 1 
and at p = for n > 2. 

The different behavior of the exchange energy in the 
various cases leads to dissimilar results. 

For n = 0, the fact that the minimum of £ Xj q occurs at 
p = 1 leads to an enhancement of p m in i-e. to the familiar 
enhancement of the spin polarization. The opposite is 
true for n > 2 since in these cases the minimum of £ x ^ 
occurs at p = 0. 

For completeness we remark that in the limiting case of 
very large winding number n the result can be obtained 
simply neglecting the cos n8 contribution in Eq. (fT7)l . For 
comparison the corresponding curve is shown as a dashed 
line in Figure [T] Notice that in this case the magnitude 
of the exchange energy is minimum for p = 1, the value 
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FIG. 1: Plot of the exchange energy per particle £ X n(r 3 ,p) 
(in Ry units) as a function of p. Here r s = 1. The different 
values of n are noted, and the limiting curve n = 00 is also 
displayed (dashed line). 



being given by: 



lim £ x , n (r s , 1) = - 



37rr s 



(18) 



The situation for the n — 1 case is more complex al- 
though in all cases the exchange only leads to a very 



4 




0.2 0.4 0.6 0.8 



y 

FIG. 2: Plot of the fractional generalized polarization p m in 
as function of the parameter 7 for different values of r s . Here 
n = 2. The increase of the depopulation value ■yd with r s is 
manifest. 



not small must be obtained numerically. As an example 
we show in Figure [21 the value of p m i n as function of the 
dimensionless coupling strength 7 for different values of 
the density parameter r s in the case of n = 2. We notice 
that for low densities the depopulation value 7^2, for which 
Pmin = 1, is considerably increased by exchange. This 
is in stark contrast with the familiar n = case. In 
particular, for the n = 2 case of Figure we have: 

Id = - + ( — — ] r s , (22) 

' 2 \12 9nJ V ' 

where the second term represents the increase of the de- 
population field due to exchange effects. In this particu- 
lar case the correction is linear in the density parameter 



small deviation from p. 



(0) 



Specifically p m in is slightly 
enhanced for p^ in < p* and slightly diminished for 



Pmin > Pi ' being unrenormalized for p^ in 



Pi- 



A similar argument leads one to conclude that the 
critical value 7^ for which the upper spin band empties 
(at fixed r s ) decreases from its non interacting value for 
n = (g = 1) while it does increase in the other cases. 

Studying the limit of small p is of particular interest 
since it corresponds to a determination of the general- 
ized susceptibility. In this case a direct inspection of the 
integral of Eq. (fTT)) leads to the asymptotic formula: 



£ ( r „)~_^_^i„2 

^x,ny stP) — „ P 

6itr s r s 
where we have defined the quantity 

r 1 



7T ' 1 — 2m 

m—0 



The resulting value for p m in is then given by: 

9 



Pn 



1 — C n r s 



(19) 



(20) 



(21) 



Eq. (|2"Tj) simply expresses the fact that in this limit the 
effect of the interactions is to renormalize the non inter- 
acting result p^ in — g via the denominator (1 — C n r s ) _1 . 
Interestingly the latter corresponds to an enhancement 
only for n = 0. 

In particular for n — we recover the well known 
Hartree-Fock differential instability occurring at r s = 
On the other hand, for n = 1 we have C\ = 
the whole renormalization effect being solely associated 
with correlations effect. Finally for n > 2, C n is (ever 
increasingly) negative leading to a perhaps iconoclastic 
exchange driven quenching of the generalized polariza- 
tion. 

Expression (|2"Tj) is valid in the limit of 7 — » or, for 
n > 1, when r s is large. The generic case when p m i n is 



IV. APPLICATION TO SPIN POLARIZED 
HOLE SYSTEMS 

Spin polarization experiments have recently been per- 
formed on GaAs two dimensional hole systems with 
growth direction along [113] and [100]. 6 In these stud- 
ies, the magnitude of the in-plane depopulation field Bd, 
when only one band is occupied, is surmised from the 
measured longitudinal magneto-resistance. While a small 
suppression of Bd with respect to its non interacting value 
B^ is observed in the [113] case, basically no suppression 
is observed for the [100] growth direction. This must be 
contrasted with the fact that, as discussed in Ref. @, in 
the absence of spin-orbit interaction and at the densities 
under consideration, the expected ratio for electrons is of 
order B°/B d - 10. 

The experiments of Ref. @ are carried out at some- 
what low densities (r s ~ 10—15) a regime in which 
a simple Hartree-Fock treatment, as well as somewhat 
more sophisticated analytical treatments designed to ap- 
proximately include correlation effects, generally fail to 
provide reliable quantitative results. On the other hand, 
Monte Carlo analysis of the system being to date non ex- 
istent, it is reasonable to expect that many of the qual- 
itative features established by a study of the exchange 
energy will prove sufficiently robust to justify such a pre- 
liminary discussion. 

In these quantum wells the introduction of a mag- 
netic field induces in the effective hamiltonian of two- 
dimensional holes the spin dependent terms given in 
Eqs. CO}, |0}, and (|A~7| of the Appendix. All these 
terms can be reduced to some of the generic forms con- 
templated by our model hamiltonian Q. The corre- 
sponding coupling strength 7 for each of these terms can 
then be extracted, as for instance explicitly done in (|A.5|) . 
(|A.6[) . and (|A.8[) . One can then make use of the defini- 
tions (JTUJ) and (fTTj) . alongside suitable numerical param- 
eters, to determine the relevant dimensionless coupling 
strength. 

For our numerical estimates we use values appropriate 
for the case of Ref. [g. In particular: W — 200 A for the 



5 



width of the quantum well, m ~ 0.2 Too for the effective 
mass, e = 12.4 for the background dielectric constant, 
and n = 3 10 10 cm -2 for the hole density (corresponding 
to r s ~ 10). 

We first consider the [113] growth direction. In this 
case the various physical properties are anisotropic in the 
plane, the principal directions being given by x = [110] 
and y = [332]. For the density under consideration the 
depopulation field is approximately given by Bd ~ 10 T 
along x and Bd ~ 5 T along y. This gives the following 
results for various dimensionless couplings: For Eq. (|A.3|) 
we have 

[^1=0.19 and 1^1=0.09; (23) 
For Eq. (fO)l : 

|fl&|=0.23 and \g v 03 \ = 0.03 ; (24) 
Finally for Eq. 

|^|=0.14 and [^1=0.16. (25) 

We recall here that (|A.3|) is a term of type n = 2, while 
(|A.4|) and (|A.7|) are both of type n = 0. We notice that 
the quadratic spin-orbit, although in general smaller, has 
a strength which is comparable to that of the terms of 
the Zeeman type. 

For the [100] growth direction, the Zeeman term (|A. 7|) 
is vanishing, and the depopulation magnetic field is ap- 
proximately Bd — 10 T. Therefore, we obtain that the 
quadratic spin orbit and the Zeeman term cubic in B 
have comparable strength: |<?03| — |<?2i | — 0.2. 

These estimates suggest that the apparent quenching 
of the many-body enhancement of the spin susceptibility 
is due to the presence of a large n — 2 spin-orbit coupling. 
Moreover the presence of a sizable n — 3 Rashba spin- 
orbit term would also result in a reduction of the gen- 
eralized polarization. This is consistent with the experi- 
mental finding that the suppression is most noticeable for 
[100] quantum wells, for which the quadratic spin-orbit 
is comparatively stronger. 

From our theory one can moreover surmise that a 
larger many-body enhancement of the susceptibility is 
expected in the limit of very low densities when, due to 
the oc k 2 and cx k 3 dependence, the n — 2 and n = 3 
terms become less relevant. The enhancement should be 
also more noticeable in the limit of a very narrow well. 
In the particular case of the [113] growth direction, this 
happens because the linear Zeeman term (|A.7[) can in 
principle become dominant. For the [100] growth direc- 
tion on the other hand, the term (|A.7[) vanishes while 
the cubic Zeeman term (IA.4I) is largest. This in spite of 
the ~ W 4 proportionality of the coupling strength. The 
reason is that the depopulation field is large in a very 
narrow well. This obtains since B d ~ 1/W* while the 
magnitude of the quadratic spin-orbit (|A.3|) behaves like 
~ Wz thereby losing its relevance. 



It should be kept in mind that, being based on a 
perturbative treatment and not taking in account or- 
bital effects, these conclusions should be taken at best 
as qualitative. 16 



V. DISCUSSION AND CONCLUSIONS 

The main conclusion of our analysis is that in the pres- 
ence of quadratic and cubic spin-orbit interactions the 
magnitude of the exchange energy decreases with increas- 
ing generalized polarization. This results in a quenched 
value of p and in a corresponding increase of the value 
of the depopulation coupling jd- By the same token the 
corresponding generalized susceptibility is also quenched. 
This interesting phenomenon stems from the in plane ro- 
tation of the spin quantization axis induced by the spin- 
orbit coupling proper (n — 1, 2, 3) and from the universal 
structure of Eq. (|15[) which only depends on the momen- 
tum space occupation (i.e. the generalized polarization 
p), and the spin orientations Sk- As one can verify by 
making use of the very same equation, for a given value 
of p the only difference between the various spin-orbit in- 
teractions stems from the different form acquired by Sk- 
Since for larger n the spins are "less parallel" it is clear 
that the magnitude of the exchange energy will decrease 
for larger n. 

Being based on a study of the exchange energy only, 
our theory can be expected to be strictly valid in the high 
density limit. The reason is that it is in this regime that 
the exchange energy represents the first interaction cor- 
rection to the non interacting result, correlation effects 
becoming comparatively smaller as r s decreases. On the 
other hand, at not too low densities, the physics of the 
exchange is still expected to give qualitatively reason- 
able results. This conclusion appears to be corroborated 
by the apparent observed reduction of the many-body 
enhancement of the spin susceptibility in dilute hole sys- 
tems in which, beside the familiar Zeeman term, the mag- 
netic field also induces large quadratic spin-orbit interac- 
tions. 

It is not difficult to prove that, when the spin-orbit 
coupling terms described by Eq. (g]) are present in isola- 
tion, the many-body states, as parameterized by r s and 
p, are also self-consistent solutions of the Hartree-Fock 
equations^ On the other hand this is not the case when 
multiple concomitant terms are present. In such situa- 
tions the circular symmetry is broken, and the interacting 
problem is considerably more complicated. Such situa- 
tions must be treated case by case. In particular the 
spin quantization directions Sk must be determined self- 
consistcntly. 17 This problem has nontrivial solutions in 
the case of broken symmetry states, like for instance in 
the case of ferromagnetic phasesi^ 

As a final remark we stress that since the generalized 
fractional polarization does not correspond directly to an 
actual magnetization, strictly speaking, one cannot draw 
direct conclusions about the enhancement of the spin- 
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spin response from measurements of the depopulation 
field B^. In general, the spin susceptibility is enhanced 
by the exchange, in a way similar to the usual case with- 
out spin-orbit X The bare spin-spin susceptibility involves 
the response to a pure n = perturbation which, from 
an experimental point of view, is not straightforward to 
realize for the case of hole systems of Ref. |g. In fact, as 
we have argued, the external magnetic field induces also 
a change of the n = 2 spin-orbit coupling. 

A detailed study of the linear spin-spin response and of 
the phase diagram, within the framework of the Hartree- 
Fock approximation and in the presence of spin-orbit cou- 
pling, will be the subject of future publications. Inter- 
estingly, the transition to a ferromagnetic state occurs at 
densities that are in general larger than in the absence 
of spin-orbit interaction^^ 



APPENDIX 

We obtain and discuss here spin dependent contribu- 
tions to the effective hamiltonian appropriate to the high- 
est two-dimensional heavy hole subband in the presence 
of an in-plane magnetic field B = B x x + B y y for a typical 
III-V semiconductor quantum well. The specific numeri- 
cal value of the parameters will chosen to be appropriate 
to the case of GaAs. 

We begin by approximately describing the motion of 
the holes in the bulk through the standard Luttinger 
hamiltonian^ which in spherical approximation takes 
the form: 



H h = — 



1 



2mn 



(71 + f() P 2 - 2 7 (J • P) 5 



(A.l) 



where jj are 4x4 spin-3/2 matrices and, for the mo- 
ment, cubic corrections have been neglected. For GaAs 
7i = 6.85 and 7 = (72 + 73)/2 = 2.5. Within this con- 
text the effect of the magnetic field can be described by 
introducing the Zceman hamiltonian Hz — — 2k hb B • J, 
where for GaAs n = 1.2, a convenient choice of the vector 
potential being provided by A = zB y x — zB x y. 

As a specific model case, we consider here the confine- 
ment associated with an infinite rectangular well of width 
W . The corresponding effective hamiltonian for holes in 
the highest two-dimensional subband can be written as: 



H =S (P) + SH 21 +5H I 



IKS 



(A.2) 



where £o(p) is the subband energy dispersion for B = 0, 
and SH21, SHqs are spin dependent terms associated with 
the external magnetic field with sub indices indicating the 
value of the integer n and the power of their dependence 
on the magnetic field. Their explicit form can be found 
by making use of perturbation theory 8 in B and in the 
wave vector k. For the first term we obtain: 



SH 2 



a/j B W 2 B + p 2 + &- + B 



(A.3) 



where B± = B x ± iB y and the explicit form of the nu- 
merical coefficient is given by a = j^s^f^^Tj^ ~ if ■ F° r 
the second term wc find: 



SHn 



bfl% 



m W 2 



B\a- +B 3 _a+ 



(A.4) 



where b 



_ k(tt -6) 



277" 



— -. — . We should remark that 

2 8(271+57) 

the present results do differ from the ones one would infer 
from the corresponding formulas appearing in Ref.[8|. For 
GaAs we have a ~ —0.2 and b ~ 1.5. 

With a suitable spin rotation, these contributions can 
be both transformed to the form of generic spin-orbit de- 
fined in our model hamiltonian ([4]). In particular, SH21 
and 5Hq3 are of the type n — 2 and n — respectively. 
Because of the isotropy implied by the spherical approxi- 
mation, the value of the corresponding coupling strength 
7 extracted by comparison to ((4]) only depends on the 
magnitude of the magnetic field B and is immediately 
found to be: 



721 



and 



703 = bn 



a^ B W 2 B 
ir 2 h 2 



, 1 W 2 m 



Tr 2 h 2 



B 



(A.5) 



(A.6) 



for 5H21 and 5i?o3 respectively. 

For a given growth direction, these results can be ex- 
tended beyond the spherical approximation to include 
the appropriate cubic anisotropy as for instance done 
in Ref. 8. Following this procedure one then obtains in 
(|A.2jl an additional anisotropic linear Zeeman term. By 
choosing coordinates along the principal axes, this can 
be generally expressed as: 

SH Q1 = ^-{g x B x a x + g y B y a y ) . (A. 7) 

As it turns out, this term vanishes (even beyond the 
spherical approximation) for the high symmetry growth 
directions [100] and [111]. However, in the case of the 
[113] growth direction, (|A.7j) is non vanishing and the 
principal axes x and y are along the [110] and [332] re- 
spectively. A perturbative estimate of the suitable Landc 
g-factors for the case of an infinite rectangular well gives 8 - 
g x ~ —0.17 and g y ~ 0.41. The coupling strength 7, as 
defined in of this n = term depends on the direction 
of the magnetic field, and is given by: 



Tor = 9 -^B 



(A. 



for the particular case of an external field of magnitude 
B along one of the two principal axes (i = x or y) . 

Finally, we mention that including a transverse electric 
field one can develop a perturbation theory in £ z and k 
to obtain a term corresponding to the n = 3 spin orbit 
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coupling appearing in Jl} . In this case the explicit form This again differs from the corresponding expression sur- 
of the coefficient reads: mised from Ref. [8|- 

512 6 ^ * (A.9) 



9?T7r 6 (71 + 2 7) (371 - 10 7) 



1 G.-H. Chen and M. E. Raikh, Phys. Rev. B 60, 4826 
(1999). 

2 D. S. Saraga and D. Loss, Phys. Rev. B 72, 195319 (2005). 

3 G. F. Giuliani and S. Chesi, in Proceedings of "Highlights 
in the Quantum Theory of Condensed Matter" (Edizioni 
della Normale, Pisa, 2005). 

4 S. Chesi and G. F. Giuliani, unpublished. 

5 Spatially inhomogeneous solutions have been investigated 
in G. Simion and G. F. Giuliani, unpublished. 

6 R. Winkler, E. Tutuc, S. J. Papadakis, S. Melinte, 
M. Shayegan, D. Wasserman, and S. A. Lyon, Phys. Rev. 
B 72, 195321 (2005). 

7 G. F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Liquid (Cambridge University Press, Cambridge, 
2005). 

8 R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, Berlin, 
2003). 

9 In this expression standard terms associated with the neu- 
tralizing background have been omitted for simplicity.—. 

Y. A. Bychkov and E. Rashba, JETP Lett. 39, 78 (1984). 

1 Y. A. Bychkov and E. Rashba, J. Phys. C 17, 6039 (1984). 



G. Dresselhaus, Phys. Rev. 100, 580 (1955). 
This is the case when the heavy hole-light hole subband 
splitting is independent of the field. See however the dis- 
cussion in B. Habib, E. Tutuc, S. Melinte, M. Shayegan, 
D. Wasserman, S. A. Lyon and R. Winkler, Appl. Phys. 
Lett. 85, 3151 (2005). 

The form of the spin-orbit coupling we consider is obviously 
not valid at large k. 

We will assume that for the n = 1 Rashba case both chiral 
bands are occupied. When this is not the case (as it can 
in principle occur at very low densities or large spin-orbit 
coupling) one can tackle the problem by introducing the 
device of generalized chirality as done in Ref. 0. 
As pointed out in Ref.@ a potentially relevant phenomenon 
we have neglected is that of the effective mass enhancement 
associated with the magnetic field. 

This also happens in the presence of terms equivalent to 
the same n by different spin rotations, as for example in 
the simultaneous presence of Rashba and Dresselhaus spin- 
orbit. 

J. M. Luttinger, Phys. Rev. 102 (1956). 



